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a series of sigma models with torsion are analysed which generate their mass dynam- 
ically but whose ultra-violet fixed points are non-trivial conformal field theories - in fact 
SU(2) WZW models at level k. In contrast to the more familiar situation of asymptotically 
free theories in which the fixed points are trivial, the sigma models considered here may 
be termed "asymptotically CFT" . These theories have previously been conjectured to be 
quantum integrable; this is confirmed by postulating a factorizable S-matrix to describe 
their infra-red behaviour and then carrying out a stringent test of this proposal. The 
test involves coupling the theory to a conserved charge and evaluating the response of the 
free-energy both in perturbation theory to one loop and directly from the S-matrix via the 
Thermodynamic Bethe Ansatz with a chemical potential at zero temperature. Comparison 
of these results provides convincing evidence in favour of the proposed S-matrix; it also 
yields the universal coefficients of the beta-function and allows for an evaluation of the 
mass gap (the ratio of the physical mass to the A-parameter) to leading order in 1/k. 
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1. Introduction 

A theory like QCD witli massless quarks in four dimensions has no explicit mass 
parameters in its classical Lagrangian; instead a mass scale A is generated dynamically at 
the quantum level. The quantity A sets the scale of low-energy physics so that the masses 
of all states in the theory, glue-balls, protons, etc., are simply numbers times A. These 
numbers are notoriously difficult to extract in QCD, either on the lattice or analytically. 
At energies much greater than A, on the other hand, the theory is asymptotically free and 
perturbation theory can be used reliably. In the language of the renormalization group 
(RG), QCD is described by a trajectory emanating from a fixed point which corresponds 
to a free theory of gluons and quarks, the direction of the trajectory being determined by 
an operator which is marginally relevant, by which we mean that it is marginal but not 
truly marginal. The fact that the operator is marginal means that no explicit mass scale is 
introduced at the fixed point itself, whilst the fact that it is not truly marginal means that 
conformal invariance is broken by the dynamical generation of the scale A as one moves 
away from the fixed point. The RG trajectory is specified by a running coupling constant 
e{h) which depends upon the mass scale scale h being probed and which in the ultra-violet 
regime (large h) behaves like 

_L.,„:,,VA).|lnln(VA).e,(lH|m). (,,) 

which in fact serves to define A precisely. In the above /3o and fii are universal numbers 
which appear as the first two coefficients of the beta-function in perturbation theory. 

A more general situation can be envisaged for a theory with dynamical mass gen- 
eration, namely the ultra-violet fixed point of the theory, while necessarily conformally 
invariant, need not be free. The purpose of this paper is to analyse such a situation in two 
dimensions in which the ultra-violet fixed point is a non-trivial Conformal Field Theory 
(CFT) - in fact a WZW model. The direction of the RG trajectory is again determined by 
some marginally relevant operator and we say that the theory is "asymptotically CFT" . 

In the case of QCD the main difficulty is the absence of non-perturbative calculational 
techniques which can be applied in the low-energy regime. In two dimensions, however, 
there is a rich class of asymptotically- free theories which are integrable: the O(A^) sigma 
models; the principal chiral models; and the Gross-Neveu models. In these theories the 
existence of higher spin conserved charges means that the S-matrix factorizes, a property 
which allows in some cases for its complete determination (see [1], [2] and [3] respectively) 
yielding an exact description of the low-energy physics. These integrable theories are 
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therefore particularly interesting from a theoretical point of view since they provide an 
arena in which one can attempt to understand the connection between the infra-red and 
ultra-violet regimes. Such a connection is also important in order to confirm the S-matrices 
written down for these models. This is because the S-matrices must, in the first instance, 
be regarded as conjectures which should be tested; in particular the question of CDD 
ambiguities must be resolved.^ 

In a series of papers ([5,6] for the O(A^) sigma model, [7] for the SU(A^) principal chiral 
model, [8] for the SO(A^) and Sp(A^) principal chiral models and [9] for the O(A^) Gross- 
Neveu model) various authors have used a technique relying on integrability to relate the 
infra-red and ultra-violet physics of families of integrable models, building on the original 
work of [10,11]. The idea is to compute a particular physical quantity - the free-energy in 
the presence of a coupling to a conserved charge - in two ways: firstly from the S-matrix 
using a technique known as the Thermodynamic Bethe Ansatz (TBA) and secondly from 
the lagrangian via perturbation theory. For the cases mentioned above the results of the 
two calculations are found to be in perfect agreement in the ultra-violet regime, thus 
resolving the problem of CDD ambiguities and, as a bonus, yielding an exact expression 
for the mass gap (the ratio of the physical mass to the A-parameter) . As well as providing 
a very stringent test of the form the S-matrix, knowing the mass-gap ratios is interesting 
in its own right as they can be compared directly with the results of lattice simulations. 

In this paper we analyse a class of theories which are asymptotically WZW models 
based on the group SU(2) at level k (see e.g. [12,13]). In accordance with the general 
situation described above, the model will correspond to an RG trajectory defined by the 
marginally relevant operator Tr(j£,J^) in the WZW theory, where Jl and Jfi are the 
usual left/right conserved Kac-Moody currents. This family of examples fits into the 
general scheme of "massive current algebras" set out in [14]. It is crucial that the theories 
we consider can also be described explicitly at the lagrangian level: they are in fact sigma 
models with a Wess-Zumino (WZ) term defined on the group manifold SU(2). This family 
of lagrangians was first written down by Balog et al in [15] who argued further that the 
resulting theories should be quantum integrable. What was not so clear in their work was 
whether the models would lie in the class of massive current algebras at the quantum level. 
Our strategy for showing that these models do lie in that class, and in particular that they 
are quantum integrable, is to use the exact S-matrices that has been proposed by Ahn et 
al [16] to describe perturbations of WZW models. We shall then use the ideas of [5-9] to 
carry out a highly non-trivial consistency check between the lagrangian formulation of [15] 
and the S-matrix written down in [16] in the manner we have already outlined above. As 
a by-product we will extract an expression for the mass gap valid to leading order in 1/k. 



^ In the case of the SU(A^) chiral Gross-Neveu model there is a derivation of the S-inatrix from 
first principles via the Bethe Ansatz [4] . 



The paper is organized as follows. In section 2 we discuss the lagrangian for the model, 
its current algebra, and its renormalization to one-loop. In section 3 we write down the 
S-matrix conjectured to describe the quantum scattering and in section 4 this is used in 
conjunction with TBA techniques to calculate the response of the free-energy to an external 
field. Section 5 contains a calculation of this same quantity in perturbation theory, after 
which we compare the expressions to confirm the choice of S-matrix and extract the mass 
gap of the model. We conclude with some further remarks in section 6. 

2. The lagrangian, current algebra and one-loop renormalization 



The integrable field theories that we shall investigate are described in two-dimensional 
Minkowski space-time (with coordinates ^^ = (r, a)) by the lagrangian density [15] 
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(2.1) 



with 



^, , cos"^ w , , I X — 1 sintycos w ,^ ^, 

(3M = ^ — ^' "("') = v^rT^ — ^- (^-^^ 

X + COS 2w \ X + 1 X + cos 2w 

The fields (WjUa) parameterize the SU(2) group manifold in such a way that a general 
group element can be written g = cos w + iUaO'a sinw where the cXa's are the Pauli matrices 
and the fields n^ are constrained via n\ + n^ + n\ = 1. e and x are coupling constants 
with X > 1. 

The complicated form of Cq requires some explanation. The most important point 
is that it ensures that the resulting theory is classically integrable - in fact it ensures 
the existence of a canonical structure consisting of two commuting current algebras [15] - 
precisely the structure studied in [16]. We shall elaborate on this point below. 

The theory has an SU(2) global symmetry generated by transformations Ua ^^ fia + 
^abcQb'i^c for parameters qi,- Finite symmetry transformations are given by the adjoint 
action g \-^ hgh~^, using some h gSU(2). This is to be contrasted with the principal 
sigma-model and WZW model which both have chiral SU(2) x SU(2) global symmetries. 
Our models are invariant under just the diagonal subgroup. 

The antisymmetric term in (2.1) is an example of a Wess-Zumino (WZ) term and as 
usual its presence leads to a quantization condition on coupling constants which is essential 
in order to obtain a consistent quantum theory. In the present case this condition is 

27r 



e^(x + l)Vx^ — 1 



/c G N. (2.3) 



One way to derive this is to consider the integral of the curl or exterior derivative of the 
WZ term over an arbitrary three-sphere, as in [12], and to demand that this always be 
a multiple of 2n. Alternatively one can require that the WZ term itself, although not 
globally well-defined, is ambiguous only up to multiples of 2%. In our case we can choose 
the ranges of our coordinates to be, for example, < w < n with Ua labelling any point on 
a two-sphere, which covers SU(2) except for one point. Then we demand that the integral 
of the WZ term should be changed by 27r on sending w ^ w + ir which gives exactly the 
condition above. ^ Yet a third possibility is to appeal to the general representation theory 
of Kac-Moody algebras because, as we shall see below, the combination in (2.3) appears 
as a central term in the current algebras which are responsible for the integrability of this 
model. 

All the information concerning the model (2.1) that we shall need in the remainder 
of this paper has now been set down. However, in view of the brevity of the presentation 
in [15] (and because there appear to be a number of numerical misprints in the relevant 
equations which can only be detected after long calculations) we shall, before proceeding, 
elaborate on the current algebra structure which is responsible for the particular form of 
the Lagrangian (2.1). We shall also supply some details of the one-loop renormalizability 
of the model which were left implicit in [15]. 

The theory (2.1) is of the general form 

/:o = ^ {G,,{4>)d^,4>'d^4P + e^''B,,{4>)d^fd,4>^} . (2.4) 

where the fields (/>*(^^) describe a map from two-dimensional Minkowski space-time to 
some target manifold. Motivated by the example of WZW models, one can ask when such 
a general sigma-model exhibits a classical current algebra. We restrict attention to the 
case in which the target manifold is the group SU(2) and it is convenient for this part of 
our discussion to choose antihermitian generators normalized so that 

^a = —-0'a, [Aa, Ab] = eafecAc , (2.5) 

which corresponds to choosing the single simple root of SU(2) to have length one. We 
shall make no distinction between upper and lower SU(2) indices. A natural Ansatz for 
the light-cone components /^ = /g ± If of a current in the SU(2) Lie algebra is 

/^ = -\Lfd+<i>' , 11 = -\R-d.<i>\ (2.6) 
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■^ The precise relationship between this criterion and the previous one is quite subtle in the 
general case; see e.g. [17]. 



where L" and i?" are vielbeins for the sigma-model metric: 

L'}L'^ = RtR'^ = G,j, (2.7) 

The equations of motion following from (2.4) ensure that these currents are conserved 
Q^ja^i, _ Q j^ pg^j^ g^igQ |-,g shown by tedious calculation that, with the canonical structure 
defined by (2.4), these currents obey a classical (equal-r) Poisson bracket algebra 



{4(a), Ilia')} = e«^^ ( a4(a) + hl^a) )5{a - a') ± ^5' [a - a') ^^ ^^ 

{lUcr),lt{a')} = -be'^'^lUcr) + r{a))6{a - a'), 

with a and b constants, provided that the quantities L" and R^ satisfy certain conditions. 
To express these conditions compactly it is convenient to introduce differential forms on 
the group manifold: 

L = A„4W, R=XaRtd(j)\ (2.9) 

Then the current algebra above will hold provided 



hL + Rh = 
dL + aL^ -hh-^L^h = Q 
dR + aR^ -bhR^h~^ =0 
3H = -aTrL^ - 3bTrRL'^ = aTtR^ + 3bTtLR 



2.-1 ^ (2.10) 



where h is some group-valued function on SU(2) and 2H = dB is the field-strength 
three-form corresponding to B. We use the conventions B = ^Bijd(j)^Ad(j)^ and H = 
^Hijkd(j)^ Ad(j)^ Ad(j)'^ for the components of two-forms and three-forms respectively; as a 
result 2Hijk = diBjk + djBki + dkBij. 

The SU(2) WZW model corresponds to a special solution of the current algebra con- 
ditions above in which 

L = g-^dg, R = -dgg-\ h = g, a = l, 6 = 0, (2.11) 

and in this case (2.8) clearly collapses to two commuting Kac-Moody algebras. The action 
written in (2.1) corresponds to a slightly more complicated solution of (2.10) which can be 
motivated as follows. First consider the Ansatz 

L = c(/i-M/i) + A„n„/(/9)dp, R = c{-dhh-^)-\anaf{p)dp, (2.12) 

where h = exp(— Aafiap), c is some constant, f{p) is a function to be determined, and 
the variable p{w) is itself some function of our SU(2) coordinate w which we will fix in 



a convenient way at the end, after finding a solution. Tfiis seems on the face of it to be 
a rather redundant procedure, but it turns out to simphfy some technical aspects of the 
discussion. The Ansatz above is clearly a straightforward modification of the WZW case, 
and it is chosen so as to satisfy the first equation in (2.10) automatically. It is not difficult 
to check that the remaining conditions in (2.10) hold if 

1 -1-1 

a = 2x + l, 6 = -l, c= f{p)= ^/°'^~ ■ (2.13) 

2(x + l) x + lx + cosp 

This solution can also be expressed in the form 

L= Xaiadua - Peabcmdric + 'jUadp) 
R = -Xaiadua + Peabcmdnc + 7^adp), 

where the functions ct, f3 and 7 are given by 



(2.14) 



sitip 1 - cos p 1 f^.^. 

" 2(0: + 1)' ^ 2(x + l)' ^ 2(x + cosp)' ^ ^ 

The sigma-model metric and WZ term are now determined as functions of p and Ua by 
the equations (2.7) and (2.10) respectively. 

The final step is to relate p to w, which can be done in such a way that the expression 
for the WZ term can be written in closed form. This is achieved by choosing 



a p X — 1 /„ i^N 

— = — cot- = -i/ tanw. (2. Id) 

/3 2 \ x + 1 ^ ' 

Using this one can deduce the expressions a{w) and I3{w) given in (2.2) and, after some 
effort, one then recovers (2.1). 

The current algebra corresponding to the solution (2.13) above was first considered 
by Rajeev [18], who showed that it could be decomposed into two commuting Kac- Moody 
algebras. The combinations which achieve this are 

obeying 

{Jl{a\ Jiia') } = e-'^JUcr)6ia - a') ± \ r^-r ^'(^ " ^') 

2e^{x + l)Vx^ -1 (2.18) 

{JS(a),4(a')} = 0. 

The ± signs occur in the central terms because these are classical Kac-Moody algebras, 
and the quantization condition (2.3) can now be recovered by comparison with some stan- 
dard reference {e.g. equation (2.3.14) of [13]). Unlike the WZW case, the components of 
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these Kac-Moody currents are not chirally conserved (although the original current /" is 
conserved by construction). Notice, however, that the WZW case can be recovered by 
taking the limit x — > oo, /c fixed, provided we rescale the fields appropriately. 

Since the theory (2.1) is a generalized sigma model (a sigma model with a WZ term) 
its renormalization group fiow can be analysed using the background field method (see for 
example [19]). We can simply quote the well-known results for the way that the metric 
and WZ term in (2.4) run with the renormalization scale to one-loop, but we must then 
ensure that these equations are indeed consistent with the specific Ansatz of (2.1). In 
our discussion of the current algbera, it was convenient to keep the coupling constant 
dependence explicit, but to apply the general renormalization results of sigma- models it is 
better to absorb the coupling constant e into our definitions of the metric and WZ term 
by defining gij = Gij/e^, bij = Bij/e^ and hijk = Hijk/e^. The coefficients of the beta- 
function are calculated in terms of the generalized curvature Rijki corresponding to the 
connection 

Tj'/c = Tj^fc + h'^jk (2-19) 

which involves the usual Christoffel connection Tj'^k (constructed from the metric gij) 
modified by a torsion term. To one loop one finds [19] that under the renormalization 
group transformation of the subtraction point /x the metric and anti-symmetric field satisfy 

where Rij = R\jk is the generalized Ricci tensor. 

We now apply these formulae to the theory (2.1). First of all, we define the coordinates 
9 and t/j via Ua = (cos 9, sin 9 cos t/j, sin 9 sin t(j) . In these coordinates the metric has non-zero 
components 

9ww = 2^2 TT' 999 = ^27 — r~iT' 9ip^ = ^2? — r~TT^^'^ ^' (2-21) 

and the anti-symmetric field has non-zero components 

1 /7r 



W ~ ~bi^6 — 



(x + l) 



V^^"^ V2 



w ] — a 



sm9. (2.22) 



From these we find that the non-zero components of the generalized Ricci tensor are 



f> 2 ^ , x + l 



x^ — 1 \ X ~ 1 

2x I 

^ee = -^^^ + 2^ol\/x'^-\, (2.23) 



R^^ = Rbb sin^ 9, 



Re^ = -R^9 = 2(3(3Wx^ - 1 sin ^ 



where the prime denotes a derivative with respect to w at constant x. Using the expression 
for the Ricci tensor in (2.23) in the equations (2.20) shows that under renormahzation 
group flow the form of the lagrangian is preserved up to a renormahzation of the couphng 
constants e and x: 

//^ = 7^(1 - 2a:)e3 + 0{e'), ^^ = -{x' - l)e' + 0{e'), (2.24) 

and a diffeomorphism of the fleld w given by 

dw 1, 9 .coswsinw n ,^, a. , . 

li— = — Or' - 1^ — e' + 0{e^). 2.25 

on TT x + COS 2m; 

These one-loop results (2.24) agree with the analysis of [15]. Notice that to this order 
k deflned in (2.3) is constant under renormahzation group flow as we expect. In the ultra- 
violet, // ^ oo, e runs to zero and x runs to inflnity. In this limit one can easily show 
that 

Co = Cwzw + ^Tr (JlJr) + Oil/x^), (2.26) 

where £vvzw is the usual SU(2) WZW lagrangian at level k, and Jl = g~^d^g and 
Jfi = —{d-g)g~^ are its left and right conserved Kac-Moody currents. This expression 
justifles our earlier statement that the theories (2.1) are SU(2) WZW models perturbed 
by the operator Tr( J^J/j). It is also easy to see that this perturbation breaks the chiral 
SU(2) X SU(2) symmetry of £wzw to the diagonal, or adjoint, SU(2) subgroup mentioned 
above. 

Assuming that k is indeed constant we can eliminate x from (2.24) to get the flow 
equation just involving e. Later we shall be interested in this flow equation for large but 
finite values of k. In this regime we deduce from (2.24) that 

^^ = -/3oe'-/3ie3-0(e4), (2.27) 



where 



(3o = \\ (1 + 0{l/k)) , /3i = -1 (1 + 0{l/k)) . (2.28) 

\ nk TV 



Notice that terms coming from higher loops can produce corrections of lower order in 1/k 
assuming that the coefficient of e^ in nide/dn) is polynomial in x. It is important to 
remember that the expressions (2.28) are universal. 

We have now described in some detail the lagrangian field theory we wish to study, 
and in the next section we conjecture an S-matrix to describe the scattering of states in 
this theory. We will subsequently undertake a non-trivial check of the form of the S-matrix 
by using the ideas of [5-9]. To do this we need to couple the theory to a conserved charge. 



The idea is to modify the hamiltonian H -^ H — hQ, where Q is a conserved charge 
corresponding to some generator of the SU(2) symmetry of the model. In the Minkowski 
space lagrangian picture this is achieved simply by replacing the derivative of n^ in the 
time-direction by the "covariant derivative" : 

doUa -^ doUa + 2heabcqbnc, (2.29) 

where the Qa are a set of parameters which later we take to be g = (1, 0, 0) without loss of 
generality (due to the SU(2) symmetry). We will then compute the response of the free- 
energy per unit volume 6f{h) = f{h) — /(O) in the ultra-violet regime in two ways: using 
the S-matrix along with thermodynamic Bethe Ansatz techniques and using conventional 
perturbation theory. 

3. The S-matrix 



Consider, for a moment, the most general way of associating S-matrices to the Lie 
algebra SU(A^). The particles form multiplets associated to the fundamental, or completely 
anti-symmetric, representations of the algebra, and each particle carries in general say m 
copies of the quantum numbers of that representation. The general two-body S-matrix 
element - from which all the others may be deduced by factorization - has the block form 
[20,21] 

StM,-M^(^) = ^"'(^)^(t)(^) ® ^(t)(^) ® ■ ■ ■ ® ^(t^)(^)' (3-1) 

where factor Sfj^ J9) acts between the j^^ copies of the fundamental representations a and 
b (a, b = 1, 2, . . . , A^ — 1) and the /cj's are parameters or coupling constants. The part 
X"'^(9) is a scalar factor which ensures that the overall S-matrix has the right analytic 
structure. Each block S^j^JO) is invariant under the action of the quantum loop-group 
Uq{SlJ{N) (g) C[6', 0~^) where the deformation parameter is q = — exp(— ^/(A + k)). In 
the limit k = oo the quantum loop-group reduces to the ordinary loop-group and the block 
'^?oo)(^) is invariant under the action of the group SU(A) itself. When k is a, natural 
number the blocks are proportional to RSOS solutions of the Yang-Baxter equation and 

^(t)(^) = l- 

The S-matrix that was proposed in [16] to describe the perturbation of the WZW 

model of level k is S^j^ ^{9). For the case of SU(2) there is only one particle and X^^{9) = 

— 1. It is worth pointing out that this general form subsumes the well-known S-matrices 

of the SU(A) chiral Gross-Neveu model, given by S^^^{e) = S^^^^^{e), and the SU(A) 

principal chiral model, given by 5"?^ oo')(^)- R-emarkably, this implies that the model (2.1) 

is equivalent at the quantum level to the SU(2) chiral Gross-Neveu and principal chiral 



models, for k = 1 and k = oo, respectively. We shall make a comment about these 
equivalences at the end of the paper. 

We now write down the S-matrix that is proposed to describe the theory (2.1). As is 
conventional, we take the kinematic variable to be the rapidity difference 9 of the incoming 
particles.'^ The S-matrix describes one massive particle with internal quantum numbers and 
for the two-body process it has the product form mentioned above [16]: 

s{e) = Ssvi2){0)(^sllf{e), (3.2) 

with k E N being identified with (2.3). The product form means that the particle carries 
two sets of quantum numbers and each factor acts on one of the sets only. The first 
factor S'su(2) {&) is the S-matrix of the SU(2) chiral Gross-Neveu model; hence the particle 
transforms in the two-dimensional representation of SU(2) and the two-body S-matrix 
elements may be written [3] 



r fi - ^) r (i + ^) 

'^S\J(2)[t') ~ 



, 9 + 2ni ,^ 



(3.3) 



^ \^ ^ 2-Ki) ^ \2 2-Ki) 

where Pt,s indicate the triplet and singlet channels. This part is equal to —Sl^^AO) as 
written above and is invariant under SU(2). 

The other factor in (3.2) describes the scattering of kink degrees-of-freedom carried by 
the particle. The particle can either be a kink or an anti-kink which interpolates between 
a set k -\- 1 vacua {l,2,...,/c-|-l} with the following selection rule: a kink can connect 
vacuum a with a + 1 and an anti-kink a with a — 1. The S-matrix is that of soliton 
scattering in the restricted sine-Gordon model [22] so the S-matrix element for the process 

Kda{ei) + Kabie2) ^ Kdcie2) + K,t{ei) is 



^kink fa ^^/^^ ^(^) /sinh(7ra/p)sinh(7rc/p)^ ^Z^'"' 



W \d c)^ ' 2^ \sm\v{nd/p)sm\v{-Kh/p) 

/sinh(7ra/») sinh(7rc/»)\ . f i'^ — d\ ^ 

^^^ • V. J/ { ■ u // + ^'^^ '^- 

\sinh[na/p)sinh{nb/p) J \ P J 



where p 




(3.4) 



RMRn{i7v-e) 

Rn{0)Rn{t7^) 



(3.5) 



We remark that the form of the S-matrix reflects the form of the lagrangian: the SU(2) 
part manifests the SU(2) symmetry of the model and the kink part describes degrees-of- 
freedom associated to the periodic field w. 



The velocity and rapidity of a particle are related by w = tanh0. 
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The free-energy from the S-matrix 



In this section we evaluate the response of the free-energy Sf{h) to the coupling with 
the charge directly from the S-matrix. The technique we use is known as the thermo- 
dynamic Bethe Ansatz [23]. In its most general form it leads to an expression for the 
free-energy in a cylindrical geometry coupled to a conserved charge which plays the role 
of a chemical potential. In our case we wish to evaluate the free-energy in the plane, i. e. 
at zero temperature. 

Consider the (one-dimensional) statistical mechanics of a gas of particles described by 
the S-matrix (3.2). Since this theory is integrable, the number of particles is conserved 
under interaction and it is meaningful to consider single particle energy levels. In a free- 
theory the energy of these levels would simply be e(^) = to cosh ^ — h, where h is the 
chemical potential, and the free-energy (per unit volume) at zero temperature would be 
that of a free one- dimensional relativistic fermion gas:"^ 



?77 

f(h) = — I dee{e)coshe, (4.1) 



where 9^ the Fermi rapidity is determined by the condition that e(±^F) = 0. In our 
case, the complications are two-fold: the theory is, after all, not free and furthermore the 
particles carry internal quantum numbers. As a result of the former e(^) now satisfies an 
integral equation involving kernels related to the S-matrix of the theory. The effect of 
the internal quantum numbers is to couple the energy e(^) to the "magnon" energy levels 
^p(^)7 P = 1, 2, . . . , /c — 1, and CqW^ g = 1, 2, . . . , oo; where the former result from the kink 
part of the S-matrix and the latter from the SU(2) part. The free-energy is then still given 
by (4.1). The equations are known as the TBA equations and they have been derived at 
finite temperature and zero chemical potential for our S-matrix in [24].^ At T = and in 
the presence of a chemical potential coupling to the charge of the SU(2) symmetry, the 



^ The fact that particles should be treated as fermions in the thermodynamic Bethe Ansatz 

results from that fact that the S-matrix satisfies 5(0) = — 1. 
^ The TBA equations for the more general S-matrices (3.1) was considered in [21]. 
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TBA equations adopt the form 

k — l oo 

e+ ie) + R* e- (9) + J^ af ) * e+ (9) + J^ «i°°^ * C (^) = ^ ^^^h 9 - h, 

p=l g=l 

fc-1 

C(^) + E 4? * ^.^(^) = 4'^ * ^"(^)' (4.2) 

9=1 

oo 



Cp"(^) + E 4r^ * C(^) = 4°^^ * ^"(^) - '^hp, 

q=l 

where star denotes the convolution / * g{9) = J d9'f{9 — 9')g{9') and f^{9) denote the 
positive/negative decomposition of f{9) = f^{9) + f~{9), i.e. f^{9) = f{9) if f{9) > or 
f{9) < 0, respectively, being otherwise zero. The kernels in (4.2) are given by 

„/^N f°^ dx ,^ , sinh (-kxIT) ,, , , 

R(9)= —cos(9x) — -— - — , \ ' ' — -exp/c7rx/2, 

^ ' Jq -K ^ ^ sinh(7r/cx/2) sinh(7rx) ^^ ' '' 

/i(fc)/m _ /"°° ^^ ta ^ ^ sinh(7^pa;/2) sinh(7r(A; - q)x/2) cosh(7rx/2) , . 

""^''Jo V'°'(^"^^ sinh(7ra;) sinh(7ra:/2) ' ^ ' 

,(fc)/m _ 1 sm{7vp/k) 



<H^)=-r 



tt/c cosh{29/k) — cos{np/k) ' 

for q> p {Apq (9) = Aqp (9)). The dependence of the free-energy (per unit volume) on the 
chemical potential is then given by 

df{h) = f{h) - /(O) = — / d9e-{9) cosh^. (4.4) 

The problem before us is to solve the coupled integral equations (4.2). Our strategy 
will implicitly assume that the solution of the equations (4.2) is unique. Given this the 
crucial observation is that ap (9) is a positive kernel; hence the solution of the TBA 
equations is Cp{0) = Cq{0) = with ^-{9) = af^ * e-{9), C~{9) = 4°°^ * e-(^) - 2hp and 

e+(^)+i?*e-(^) =mcosh^-/i. (4.5) 

The solution e(^) to (4.5) is some symmetric convex function with zeros at the Fermi 
rapidity ±9-F{h). When h < m the system is below threshold; the external field is too 
weak to excite any particle states and hence 6f{h) is zero. Beyond the threshold h = m 
the chemical potential forces the system into a state where the particles line up their spins 
with the external field. From the form of the TBA equations we see that the external field 
does not couple to the kink number and it turns out that the ground-state has total kink 
number zero. 
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Notice that the reasoning which led to identifying the solution of the TBA equations 
in terms of a particular configuration of the quantum numbers of the particles was arrived 
at from studying the full TBA equations. This is to be contrasted with the more heuristic 
arguments used in [5-9] leading to the hypothesis that only one particle-state contributed 
to the ground-state. 

We have arrived at an expression for the free-energy in terms of a single function e(^) 
which satisfies the single integral equation (4.5). This equations are of the same form, but 
with a different kernel, as those of the O(A^) sigma model [5,6], principal chiral models 
[7,8] and Gross- Neveu models [9]. 

It is not possible to solve the equation (4.5) in closed form; however, we will be 
interested in the solution only in the deep ultra-violet h ^ m for which one can develop a 
series solution using generalized Wiener- Hopf techniques [5,25] (see the appendix of [9] for a 
clear summary) . Rather than explain these techniques we simply follow the manipulations 
of [9] required to extract the series solution. 

The method starts by decomposing the Fourier transform of the kernel R{0) in the 
following way: 

sinh (7rx/2) ., , x 1 , , . 

^ /^ exp(A;7ra:/2) = — ^— ^^^, (4.6) 



sinh(7rA;x/2) sinh(7ra:) G+{x)G-{x) 

where G-{x) = G+{—x) and G±{x) are analytic in the upper/lower half-planes, respec- 
tively. So 

^ / ^ /7T7 T^l-ix/2) f. ^ .kx, , . A ,,^, 

^+(^) = ^' V(l-zfcx/2)r(l-zx) ^"P [^xh-^-H-^x)^ , (4.7) 

where 

6=^_ln2--ln-. (4.8) 

2 2 2 ^ ' 

Following [9] we now define the function a{x) = exp{2ix9F)G - (x) / G +{x) , where €{±9^) = 

0. a{x) has a cut along the positive imaginary axis and we define 7(^) in terms of the 

discontinuity: 

a{i^ + 0) - a{i^ - 0) = -2ie"2^^^7(0, (4.9) 

giving in this case 

n^j ^^ ^ ^ ^ ^ r2 (1+^/2) r(l-A;^/2)r(l -0 ^ ^^ ^ ^ ^ 

If one consults [9] then it soon becomes apparent that 7(^) for our model has the same 
functional form as a fermion model, rather than a sigma model, namely 



^(0 = 7re-'=«'"«5^rf.r- (4.11) 



n=l 
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The expansion of the free-energy is given in terms of the quantities dj . It turns out that 
df{h)/h'^ is power series in the effective couphng u = u{h) defined through 



1 

u 



klnu 



1 

z' 



(4.12) 



where 



1 

z 



In 



/i2 f2G+{0) 



m^ 



(4.13) 



Putting these expression together with the results of [9] allows us to extract the first few 
terms in the expansion of the free-energy as a function of h/m 



Sf{h) 

27r 



G+(0)2 {1 - 2diz + 2kdiz^ \ivz-2 [2di - T'{2)kdi - dj + rfs] 
-2k'^diz^ In^ z + 2k [4di - T'{3)kdi - 2dl + 2^2] ^^ In^ + 0{z^)] . 



(4.14) 



In the above, Sf{h) is an expansion in terms of the form z'^ln"' z where m > n. From 
(4.10) we find 

di = ^, rf2 = -Hn2-^lnQ)-F^r'(2). (4.15) 

So putting everything together for the models that we are considering we find that the 
free-energy extracted from the S-matrix yields the following expansion in the ultra-violet: 



Sf{h) 



K^k 



TV 



h h'^ h 

1 s^ s ins 

2 4 2 



■ 1 fi4/c k k k' 

1 in in — 

2 TT 4 2 4 



— s'^ln^'sH 

8 2 



1 fi4A; k k k 
1 In In — 

2 TT 2 2 4 



s^lns + C(s' 



(4.16) 



where s ^ = ln(/i/m). 



The free-energy from perturbation theory 



In this section, we develop the expansion of the free-energy Sf{h) in perturbation 
theory. We assume, following the discussion to one-loop in section 2 and [15], that under 
renormalization group fiow k is constant. Therefore, we may express 1/x in terms of e: 



1 



X 




2tc 4n ^ 



(5.1) 



14 



and we shall find that in the ultra-violet e runs to zero and hence x runs to infinity. The 
loop expansion parameter is e^(x + 1)^ which is expressed in terms of e as 

e^(. + l)^^f (l + yAe + ±e^ + 0(e3)). (5.2) 

This means that the contributions from higher loops can lead to terms of the same order in 
e, but their coefficients will be suppressed by higher powers of 1/k. Our result to one-loop 
will therefore be valid in the large but finite k limit. 

The Minkowski space lagrangian in the presence of the chemical potential is given by 
substituting (2.29) in (2.1) which on Wick rotating to Euclidean space becomes 

^ = ^°- e^(^ + l) (^-^i) 

+ — — —p-\p{w){n2don3-n3don2)+ —^===:(--w)~a{w) diui} . 
e{x + 1) [ y\/x^ - 1 V2 / J J 

The ground-state of the system is given by ni = ty = (modulo some discrete ambiguity 
depending on the precise way we parametrize SU(2)). We wish to calculate the change in 
the free-energy per-unit- volume 5/(/i) to one- loop, so it suffices to expand the lagrangian 
to quadratic order around the ground state with (n2,n3) = A/r^-n^(cos -i/;, sini/;). On 
suitably re-scaling the field w we find 



(5.4) 



-4/^2+ 4/^2^2 + 4/l2«;2/'^^^ +■ 

\X + \ ) 

where the ellipsis represent the interaction. 

We can simply read off the tree-level contribution to 6f{h): 

e^[x + ly TT I V 27r 

The one-loop contribution is, from the quadratic Euclidean lagrangian written above, 

where the operator M acts on (ni, w) according to 
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Here // is a mass-scale introduced to make the eigenvalues of M dimensionless. Note also 
that the field -^ is completely decoupled to this order in the loop expansion. 

The operator above is rather unconventional in nature and it proves convenient to 
evaluate its determinant using zeta-function techniques. Some basic facts concerning these 
methods, together with their application to operators of the above type, are summarized 
in the appendix. The important point for our purposes is that 

lnDetM = -CV(0), (5.8) 

where the zeta-function C,m{s) corresponding to the operator M can be represented by 

with 

Ah'^ix'^ + 1) Sh'^x Ahx _ ^^, 

The factor V denotes formally the volume of two-dimensional spacetime (strictly speaking 
this should be dealt with using some explicit infra-red regularization but these details 
are irrelevant for our purposes). Unfortunately the integral above cannot be evaluated in 
closed form, but we shall sketch below how it can be successfully expanded in powers of 
1/x to the order which we need. (Recall that the ultra-violet limit of our models is a; ^ oo 
with k fixed.) 

The strategy is to expand the cosh factor in (5.9) as a series and collect terms of a 
given power in pi so that the momentum integrals can be evaluated. The t integrals can 
then be expressed in terms of F-functions, and one obtains the result 

Cm{s) = ^^ Yl CmX-^'+-^r{s + m) , (5.11) 

Tn= — 1 

where the coefficients are given by the rather complicated expressions 

^-- 2^ (2n)!(r!)2(n-r)!22-^ ^ ' ^^'^^^ 

m = 2n-r-l \ / \ J \ / 

It is not obvious how to sum the series above completely, but we note that p in these 
expressions is 0{l/x) which will enable us to simplify things presently. 
On differentiating and setting s = we find 



Cm(0) 




(2 \ \ oo ^ 

A - ^ j In — + 5^ C^\-^V{m) \ . (5.13) 
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At this stage our result is still exact, but we must now consider its behaviour as a power 
series in p (and hence 1/x) to make further progress. The leading contribution is 0{p^); 
extracting all such terms from the expressions (5.12) for Cm with m>l gives us an infinite 
series which can be summed^ to yield 

The next contributions are 0{p^) and again the series resulting from the coefficients Cm 
with m > 1 can be summed^ to give 

The other contributions to (5.13) are O(p^) = 0{l/x'^) and we neglect them. Now we 
substitute (5.14) and (5.15) in (5.13) and expand each of A, p and rj in powers of 1/x using 
(5.10). It turns out that although (5.14) is 0{p'^) and p = 0{l/x) the particular form of 
the expressions for A and t] imply that (5.14) is, more precisely, 0{l/x'^). The final result 
is 



Cm(0) 



2 



V Ah 






+ 0{l/x^). (5.16) 



Now to obtain the required determinant in (5.6) all we need do is restore the factor 
e^(x + l)^ (which we dropped in (5.8) for simplicity) by rescaling p and use the consequence 
(5.1) of the quantization condition (2.3) to eliminate x in favour of e and k. This gives a 
final result for the change in the one-loop contribution to the free-energy per unit volume 

W».^^{f(l-.n5^)+0(e3)|. (5.17) 

The change in the free-energy dfo + 6fi which we have just calculated is a renormal- 
ization group invariant quantity, i.e. it is independent of p when the coupling constant 
runs with p. This fact allows us to determine the running of the coupling, as expressed by 

the beta-function 

de FY 



although the coefficient Pi is not determined to the order that we are working. Since 
the first two coefficients of the beta-function are universal, it is comforting that (5.18) is 
consistent with the calculation of the beta-function via the background field method (2.28). 



Er=2 M^)y^ = (1 - y)(l"(l - y) - 1) + 1 for suitable y. 

■>oo 1 ^2n— 1 ]_ , 

,n=l 2n-l^ ~ 2 ' 



^ Er=i 2;h[^^"'^ = l(ln(l + z)- ln(l - z)) for suitable z. 
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It is convenient to use the fact that the free-energy is independent of // to set // = /i and 
then, integrating (5.18), one obtains 



l/^e + /?i(^)'^^^lne + /??(^)'^'e^(ln^e + ln0+O(a, (5-19) 



where ^~^ = ln(/i/A^) is defined in terms of the A-parameter of the zeta-function regular- 
ization scheme and the expansion is in terms of the form ^"^ In" ^ with m > n. 
Hence to one-loop we deduce 






-r in^ e - 



(^ + l-ln^)^^ln^ + 0(a} 



(5.20) 



The effect of higher loops would be to introduce corrections at the same order in ^ but 
with coefficients which are suppressed by powers of 1/k. 

We now reach our main result: at the order to which we are working the two expansions 
(4.16) and (5.20) are consistent, a fact which provides a highly non-trivial check of the 
conjectured S-matrix. Furthermore, by comparing the two expressions we can extract the 
mass-gap ratio for large k: 

In ^ = -1 + ^ ln2 + - + - In - + 0{l/k). (5.21) 

A^ 2 2 4 2 4^'^ ^ ^ 

We also deduce that the second coefficient of the beta-function is, for large /c, simply 
Pi = — I/tt + 0{l/k), a result which is in perfect agreement with the second coefficient of 
the beta-function computed directly using the background field method (2.28). 



6. Conclusions 

We have investigated a series of theories that generate their mass dynamically but 
which are asymptotically (in the ultra-violet limit) non-trivial CFTs. The theories are 
in addition integrable, a property implying that their S-matrices factorize, allowing us 
to conjecture a form for these S-matrices. The calculations of the free energy which we 
carried out provide a highly non-trivial check on the form of the S-matrices since, as pointed 
out in [7], the addition of any CDD factors would drastically alter the thermodynamics 
of the system and consequently destroy the remarkable consistency between the TBA 
calculation and the perturbative result. We can conclude with some confidence therefore 
that the proposed S-matrices are correct and that the classical integrability of these models 
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extends to the quantum regime. Notice that the S-matrix has a quantum group symmetry, 
an invariance which does not seem to be manifested at the lagrangian level in any simple 
fashion. 

It is worth pointing out that the leading order behaviour —h^k/n of the free-energy 
near the ultra-violet fixed-point can easily be deduced from the knowledge that in this limit 
the theory is an SU(2) WZW model of level k. At the fixed point the chemical potential 
couples to the combination J = Jl + Jr, where Jl and Jr are the left and right currents 
of a U(l) subalgebra of the SU(2) current algebra. Following [26], the response of the 
free-energy in a finite volume V is given by 

where the expectation value is evaluated in the WZW model. The final result is propor- 
tional to the anomaly in the U(l) currents, however one must take careful account of the 
form of the operator products in a finite volume [26] : 

~ ~ k 2'Kk 

{Jl{z)Jl{w)) = -^ + -— , {Jr{z)Jr{w)) 

[z — wy V 

{Jl{z)Jr{w)) = 2nk6^^\z - w) - 
It is then straightforward to to extract [26] 





k 


+ 


27Tk 


(z- 


-wy 


V 


2nk 









(6.2) 
V' 



Sf{h) = - — , (6.3) 

IT 

in agreement with the leading order term in (4.16) and (5.20). Obviously one could extend 
this calculation away from the fixed-point by perturbation theory and hope to reproduce 
the series (5.20); a calculation which would be interesting since it would be valid not just 
in the large k regime. 

We should also emphasize a remarkable consequence of the equivalence of the la- 
grangian and S-matrix descriptions we have established in this paper: namely that the 
field theory (2.1) for /c = 1 and /c = oo is quantum equivalent to the SU(2) chiral Gross- 
Neveu model and principal chiral model, respectively (since our Ansatz for the S-matrix 
then reduces to these well-known cases). Let us consider the latter equivalence in more 
detail. In taking the limit k -^ oo, at fixed e, it is necessary to introduce the field 
r = ((7r/2) — w)/\^x'^ — 1. The lagrangian then has a well-defined limit 

1 ( 2 r^ 2 2r^ 1 

^°"2^V^''''^ + 1 + 4r2 (^^"""^ + ^ ^ ^^^ ^abce'^'nad^nbd^n.j . (6.4) 
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With a simple change of variables (j)a = rria, this may be written 

/:o = ^ (iTl^) ^^'^"^ + ^'^"'^'^ d^(l>ad^h + 2eabce^''(t)ad^,(t)bd^(t)c} ■ (6.5) 

This lagrangian is actually the non-abelian dual of the SU(2) principal chiral model [27] 
and hence it is indeed known to be quantum equivalent to it. It would be interesting to 
show, in a similar spirit, that (2.1) with k = 1 was a bosonized form of the SU(2) chiral 
Gross-Neveu model. 

Finally, it would clearly be interesting to extend the various results above to larger 
groups and also to S-matrices of the more general form (3.1). 

TJH would like to thank Michel Bauer for many interesting discussions. 



Appendix : Some basic facts about zeta-functions 



The zeta-function C,m associated with an operator M with eigenvalues A^ can be 
thought of formally as 

i 

It is not immediately obvious when this formula makes sense beyond the simplest circum- 
stances in which M is finite-dimensional with A^ > and Re(s) sufficiently large. By 
comparison with this simplest case, however, it should seem reasonable that when C,m 
exists it can be used to calculate the dimensions and determinant of M via 

dim M ^Cm (0) , Det M = -Cm (0) • ( A.2) 

To make use of this we need some means of calculating Cm (without first finding all its 
eigenvalues and using (A.l)!) and this is provided by the heat-kernel representation. The 
key is the observation that (A.l) is formally equivalent to 

1 f°° 
Cm{s) = -— dtt'-'Trg{t), (A.3) 

r(s) Jo 

where Q{t) = exp(— tM). Then we note that Q{t) = exp(— tM) can also be characterized 
by the 'heat' equation and boundary condition 

^^-Mg, ^(0) = l, (A.4) 

(where in the boundary condition 1 is the identity operator on the relevant space). Now 
we can dispense with (A.l) entirely by adopting (A.3) and (A.4) as a definition of Cm{s) 
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for any given (suitably well-behaved) operator M and we can take (A. 2) as providing 

definitions of the dimension and determinant of M. 

The case of interest for us is an operator in two Euclidean space dimensions of the 

form 

1 / —9^ + a? cdi 



A solution of (A. 4) can be found by Fourier transforming to momentum space: 

After substituting in (A. 3) we can evaluate the functional part of the trace immediately 
to obtain a factor V which is the two-dimensional space-time-volume. This leaves us still 
with a matrix trace 

rj^^^-tMip) ^ 2e-(*/A.^)(P^+A) cogh (^{t/^^){cV, + p2)V2J ^ (A.7) 

where 

A=l(a2 + 62), p=l(a2-62). (A.8) 

On combining this with (A. 3) and (A. 6) and re-scaling t we obtain the expression (5.9) in 
the text with tj = c. 
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